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Abstract 


This  paper  presents  the  methodology  for  a  system  reliability  analysis  that  can  determine  the  probability 
of  mission  success  of  the  destruction  of  buried  concrete  targets  using  deep  penetration  weapons.  The 
success  of  the  mission  depends  on  various  parameters,  such  as  the  properties  of  the  target,  the  velocity  and 
structural  strength  of  the  missile,  the  impact  angle,  etc.,  as  well  as  other  military  scenarios.  However,  this 
paper  only  considers  the  structural  design  aspects.  In  this  paper,  analytical  equations  for  the  depth  of 
penetration  and  buckling  strength  of  the  missile  are  used  to  demonstrate  an  efficient  system  reliability 
analysis  methodology.  These  equations  are  also  used  to  determine  the  sensitivity  of  mission  success  to 
various  parameters.  A  discussion  on  the  selection  of  parameters  to  maximize  the  mission  success  is  also 
presented. 

©  2005  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

Most  strategic  targets  are  buried  underground  and  their  destruction  is  important  for  the  overall 
success  of  a  war.  Destruction  of  these  targets  with  a  minimum  number  of  missiles  is  economical  to 
the  military,  more  lethal,  and  beneficial  from  a  strategic  point  of  view.  The  destruction  of  these 
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targets  is  usually  measured  through  the  damage  caused  by  the  striking  missiles.  Usually  a  full 
penetration  of  the  missile  into  the  target  is  considered  as  the  complete  destruction  of  the  target 
and  would  represent  a  successful  mission.  Investigators  in  the  past  have  proposed  various 
empirical  [1,2],  analytical  [3],  and  numerical  [4,6]  procedures  for  the  estimation  of  penetration 
depths  under  a  given  missile  impact.  These  predictions  were  deterministic,  whereas  the  actual 
penetration  problem  has  significant  uncertainties  in  its  variables.  This  shows  that  the  missile 
penetration  into  a  target  is  a  random  event,  while  the  past  work  [7,8]  shows  only  a  deterministic 
approach.  In  this  paper,  a  detailed  system  reliability  analysis  of  these  missions  is  performed  by 
considering  the  depth  of  penetration  and  the  structural  integrity  of  the  missile  as  the  criteria.  The 
various  expressions  required  to  determine  the  depth  of  penetration  in  crater  and  tunnel  regions 
were  derived  based  on  the  methodology  obtained  in  the  literature  [1,5].  These  penetration  depths 
have  been  used  in  conjunction  with  the  buckling  criteria  of  the  missile.  The  sensitivity  of  mission 
success  to  various  uncertain  parameters  is  examined  and  the  important,  controllable  parameters 
are  identified.  A  robust  configuration  for  maximizing  the  mission  success  is  presented  that  is  based 
on  these  parameters. 


1.1.  Reliability  analysis 

A  fundamental  problem  in  structural  reliability  analysis  is  the  computation  of  the  multi¬ 
dimensional  integral 

Pf=  I  fxiX)  dX,  (1) 

JQ 

where  fx(X)  denotes  the  joint  probability  density  function  (PDF)  of  the  vector  of  basic  random 
variables  and  X  —  (xi,X2, . . .  ,xn)T  represents  uncertain  quantities  such  as  loads,  geometry, 
material  properties,  and  boundary  conditions.  This  equation  is  also  called  the  convolution 
integral  because  convolution  of  the  individual  PDF  over  the  failure  domain  would  result  in  the 
joint  PDF  of  the  response.  Furthermore,  Q  is  the  failure  region  modeled  by  the  limit-state 
function  or  the  performance  function  g(X).  The  failure  region  is  defined  by  g(  X)  ^  0  and  pr  is  the 
probability  of  structural  failure. 

In  practice,  it  is  usually  difficult  to  construct  the  joint  density  function  fx(X)  and/or  to 
determine  each  individual  density  function  due  to  the  scarcity  of  statistical  data.  Even  in  the  case 
for  which  the  statistical  information  may  be  sufficient  to  determine  these  functions,  it  is  often 
impractical  to  perform  the  multi-dimensional  numerical  integration  over  the  failure  region  Q. 
These  difficulties  have  motivated  the  development  of  various  reliability  approximation  methods 
within  the  past  2  decades. 

Among  the  approximation  methods  available,  the  most  common  one  is  the  first-order  reliability 
method  (FORM)  [9].  In  this  method  the  limit-state  is  approximated  by  the  tangent  plane  at  the 
most  probable  failure  point  (MPP).  The  MPP  is  a  point  on  the  limit-state  in  the  normalized 
domain  that  has  a  maximum  probability  density.  FORM  is  most  accurate  when  the  random 
variables  are  normally  distributed  and  the  limit-state  function  is  nearly  linear  in  the  neighborhood 
of  the  MPP.  However,  if  the  failure  surface  has  a  high  nonlinearity,  which  is  the  case  in  depth  of 
penetration  situations,  the  failure  probability  estimated  using  the  safety  index  ft  by  FORM  gives 
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inaccurate  results.  In  order  to  improve  the  accuracy,  high  fidelity  methods  [10-12]  have  to  be 
applied. 

The  failure  probability  of  the  system,  however,  is  the  integration  of  the  joint  PDF  over  the 
entire  failure  region  obtained  by  the  intersection  of  all  the  limit-states,  as  opposed  to  a  single  limit- 
state  in  component  reliability  analysis.  Structural  systems  are  classified  as  series  systems  and 
parallel  systems  [13].  A  series  system  is  one  in  which  if  even  one  component  fails,  the  whole  system 
fails.  These  systems  are  also  called  weakest-link  systems.  In  the  case  of  parallel  systems,  the  system 
survives  even  if  one  or  more  of  the  components  have  failed.  The  system  fails  to  function 
satisfactorily  only  when  every  component  has  failed.  Parallel  systems  are  also  called  redundant 
systems.  In  the  present  work,  a  series  system  is  considered  in  which  the  mission  fails  if  the  any  of 
the  criteria  is  not  satisfied. 

1.2.  Solving  the  convolution  integral  in  intervals 

For  solving  the  convolution  integral,  the  function,  based  on  which  integral  is  evaluated,  must  be 
available  as  a  closed-form  and  separable  equation.  The  procedure  for  solving  the  convolution  is 
discussed  in  detail  by  Penmetsa  and  Grandhi  [12].  As  a  single  approximation  cannot  model  a 
highly  nonlinear  joint  failure  region,  several  function  approximations  are  constructed  over  the 
entire  domain  represented  by  the  uncertain  variables. 

Let  us  suppose  that  the  failure  probability  of  the  structural  system  is  governed  by  two  limit- 
state  criteria,  g i(vi,X2)^0  and  g2(Ai5*2)^0,  where  x\  and  V2  are  the  random  variables,  as  shown 
in  Fig.  1.  These  limit-states  intersect  at  a  point  where  x\  =  xq.  The  joint  failure  region  comprises 
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Fig.  1.  Convolution  in  intervals. 
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of  the  part  of  g2  where  x\  <xo  and  the  part  of  g]  where  xi  >x2.  Based  on  this  information,  the 
joint  failure  region  is  divided  into  two  different  integrals,  as  shown  in  Eq.  (2). 

r  r  r*o  r° o 

/  fx(X)  dX=  /  /  fx{X)dx2dxx 

J  J  Q  J — oc  J — oc 

pOQ  POO 

+  /  /  fxiX)  dx2dx\.  (2) 

J  X  0  J  — oo 

The  two  integrals  in  Eq.  (2)  can  be  solved  separately.  The  PDF  of  x\,  when  x\<xo,  is 
convoluted  with  the  PDF  of  X2,  based  on  g2 ,  to  evaluate  the  first  integral,  since  only  g2  represents 
the  system  failure  surface.  The  PDF  of  x\,  when  x\  >xo  is  convoluted  based  on  gx  for  the  other 
integral.  The  final  PDF  of  the  system  is  a  combination  of  these  two  integrals. 

In  a  problem  with  more  than  two  random  variables,  the  intersection  of  the  limit-states  is  not  a 
point  but  a  surface  (hyperplane),  making  it  difficult  to  split  the  PDFs,  as  shown  in  Fig.  1. 
Therefore,  the  entire  design  space  is  divided  into  regions  and  the  convolution  integral  is  solved  for 
each  region  separately.  The  convolution  integral  is  then  written  as  the  sum  of  integrals  over  all  the 
regions  as 


fx(X)dX  = 


+ 


fx(X)dX 

fx(X)dX  +  ... 


02  ' 


(3) 


As  shown  in  Fig.  2,  the  entire  design  space  is  divided  into  several  regions  such  that  the  joint 
failure  region  can  be  modeled  with  approximate  models  that  are  accurate  within  that  region.  The 
convolution  integral  is  divided  into  several  integrals  as  shown  in  Eq.  (3).  These  integrals  are  then 
solved  using  fast  fourier  transforms  (FFT),  based  on  the  approximate  model  that  is  accurate 
within  that  region.  Finally,  all  the  PDFs  of  the  various  models  are  added  to  result  in  the  PDF  of 
the  structural  system. 


Fig.  2.  Multiple  approximations  modeling  different  failure  regions. 
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1.2.1.  Numerical  example  for  system  reliability 

A  cantilever  beam,  as  shown  in  Fig.  3,  subjected  to  a  tip  load  of  444.82  N,  is  selected  to 
demonstrate  the  system  reliability  method  used  in  this  paper.  Two  failure  criteria  for  the  structure 
were  considered:  (i)  the  displacement  at  the  tip  of  the  beam  should  be  less  than  0.002  m,  as  shown 
below,  and  (ii)  the  stress  in  the  beam  should  be  less  than  34  MPa. 

Displacement  limit-state  gfX)  =  —  0.002^0 

Stress  limit-state  g2(^0  =  —  34  x  1 06  ^  0 

where  L,  b,  h  are  the  length,  width,  and  height  of  the  beam,  which  are  all  taken  as  random 
variables.  The  length  of  the  beam  is  assumed  to  be  normally  distributed  with  a  mean  value  of 
0.762  m  and  a  standard  deviation  of  0.0762  m.  The  width  and  height  of  the  beam  were  both 
assumed  to  be  normally  distributed  with  mean  values  of  0.0635  m  and  standard  deviations  of 
0.00635  m. 

As  the  limit-state  functions  are  available  as  a  closed-form  expressions  in  terms  of  the 
random  variables,  there  is  no  need  to  construct  approximations  of  each  of  the  limit-state 
functions.  Using  these  closed-form  expressions,  several  points  were  sampled  on  the  joint  failure 
region.  These  sampled  points  were  used  in  the  construction  of  multiple  response  surface  models  to 
represent  the  joint  failure  region.  The  design  space  is  divided  into  smaller  domains  such  that  the 
response  surface  model  captures  the  joint  failure  region  accurately  within  that  domain.  The 
convolution  integral  is  solved  over  the  entire  design  space  to  obtain  the  failure  probability  of  the 
system. 


Table  1 

Comparison  of  results  obtained  with  Monte  Carlo  simulation 


Methodology 

System  failure  probability 

%  Difference 

Monte  Carlo  simulation 

0.0091 

— 

RSM  +  FFT 

0.0093 

2.19 

First-order  series  bounds 

0.0089-0.0111 

-2.19-21.98 
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Table  1  illustrates  the  accuracy  of  the  methodology  in  the  prediction  of  the  failure 
probability  of  the  system.  The  error  in  the  failure  probability  estimated  by  the  proposed 
methodology  to  that  of  Monte  Carlo  simulation  is  around  2.19%.  One  million  evaluations  of  the 
limit-state  functions  were  used  in  the  estimation  of  the  failure  probability  using  the  Monte  Carlo 
simulations. 


2.  Expression  for  depth  of  penetration  and  buckling  load 

Depth  of  penetration  is  a  highly  nonlinear  expression  that  depends  on  various  parameters, 
such  as  missile  velocity,  material  properties,  angle  of  impact,  etc.  Various  assumptions  have 
been  used  to  derive  the  expression  for  depth  of  penetration  based  on  the  information  available  in 
the  literature.  It  is  assumed  that  the  impact  is  normal  and  axi-symmetric,  and  that  the  missile  does 
not  carry  a  warhead.  This  ignores  the  explosion  condition  in  determining  the  depth  of 
penetration.  Moreover,  the  missile  is  considered  as  rigid  in  order  to  determine  the  depth  of 
penetration;  since  this  is  not  the  case,  the  buckling  constraint  is  added  to  restrict  maximum  force 
on  the  missile. 

A  missile  that  has  an  initial  velocity  vo  penetrates  into  the  geomaterial  experiencing  an  axial 
resting  force  FA  at  the  nose,  due  to  the  soil  medium.  This  force  for  a  given  missile  configuration 
and  geomaterial  lay-up  (Fig.  4)  were  obtained  from  Ref.  [5]. 

Fa  =  y.na2vA  tan (6)fKpg,  (4) 

where  a  is  the  slope  of  the  linear  fit  to  the  curves  plotted  between  ar/K  and  vA  tan (0)/ yJKpg,  a  is 
the  missile  shank  radius,  vA  is  the  axial  velocity  of  the  missile  at  any  time,  6  is  an  equivalent  half 
cone  angle  equal  to  tan-1  ( 1  / -JA'V  —  l)  ,  W  is  the  caliber  radius  of  the  missile  nose,  K  is  the  bulk 
modulus  of  geomaterial,  and  pCJ  is  the  initial  mass  density  of  geomaterial. 

By  using  Newton’s  second  law,  the  equation  of  motion  of  the  missile  with  mass  m  can  be  written 
as  m(d2Zg/dt 2)  =  —Fa,  where  Zg  is  the  depth  of  penetration  at  any  time  t.  By  integrating  this 


Fig.  4.  Generic  model  of  a  deep  penetration  problem. 


1448 


R.C.  Penmetsa  /  International  Journal  of  Mechanical  Sciences  47  (2005)  1442-1454 


equation  of  motion  with  initial  conditions  we  can  obtain  the  following  expressions: 
mv  o 


c xnta 2  tan (6)^Kpg 

X  .  1  _  Exp  - anta 2  tan(0), /KfTg 


m 


d  za  _ 

va  =  =  v0  Exp 


—ointa2  tan {0)yJKpg 


m 


(5) 

(6) 


By  using  Eqs.  (5)  and  (6)  we  can  obtain  the  velocity  of  missile  (r,)  striking  the  concrete  target 
where  Zg  =  dg. 


Vs  =  Vo 


1  - 


dgCtnta 2  tan (Q)yJKpg 
mv o 


(V) 


When  a  missile  impacts  a  buried  target  with  a  velocity  va  it  creates  a  conical-shaped  crater 
region  with  a  depth  approximately  equal  to  twice  the  missile  shank  diameter,  4a,  followed  by  a 
cylinder-shaped  tunnel  region  with  a  diameter  nearly  equal  to  the  shank  diameter,  2a,  as  shown  in 
Fig.  5.  Expressions  for  these  depths  are  given  below: 


Z crater  — 


na 


1  - 


dgOtna 2  tan(0)  yJKpg  1 
mv  o  | 


m  +  47ia3Npc 


^  Sfc  +  Npcvl  1  —  (dgvina2  tan (6)^/Kpg  ) /mv o 


(8) 


4a 


◄ - ► 

2a 


Fig.  5.  Cartner  and  tunnel  formation  for  a  group  penetrating  missile. 
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Z tunnel  —  4$  4 


m 


x  In  1  + 


2na2Npc 

Npc(m2vl 


2mvodgOL7ia2  tan {6)^Kpg  +  d2gd2Ti2a4  tan2  (9)Kpg  —  4na3mSfcy 


( m2Sfc  +  4na3mNpcSfc) 


(9) 


where  dg  is  the  depth  of  target  in  the  geomaterial,  fc  is  the  compressive  strength  of  concrete,  pc  is 
the  density  of  concrete,  and  N  and  S  are  dimensionless  constants.  By  using  these  expressions,  the 
depth  of  penetration  of  the  missile  can  be  obtained  as  follows: 


Missile  depth  of  penetration  G\ 


If  % crater  4a  then  Gi  —  Z tunnel •> 
Else  Gj  —  zcraier. 


(10) 


The  second  limit-state  is  the  buckling  constraint  on  the  missile.  The  missile  is  assumed  to  be  a 
flexible  structure  and  its  critical  buckling  load  is  modeled  using  the  following  equation: 

E(na4/4)n2 


PCr  = 


(1.2  Lf 


(11) 


By  using  the  maximum  force  expression  from  Eqs.  (4)  and  (11),  we  can  easily  formulate  the 
second  limit-state  equation  to  maintain  structural  integrity. 


3.  Problem  definition 

For  the  purpose  of  performing  a  system  reliability  analysis,  a  missile  of  mass  182  kg  and  shank 
radius  0.0825  m  with  impact  velocity  41 1  m/s  has  been  selected.  These  data  are  obtained  from  Ref. 
[5].  Material  properties  of  the  geomaterial  and  concrete  are  given  in  Table  2.  The  selected 


Table  2 

Uncertain  variable  information 


Random  variable 

Distribution 

Mean 

cov 

Reference 

Mass  of  the  missile,  m 

Lognormal 

182  kg 

0.05 

[5] 

Missile  shank  radius,  a 

Normal 

0.0825  m 

0.05 

[5] 

Caliber  radius  head,  W 

Normal 

9.25 

0.05 

[5,2] 

Missile  impact  velocity,  v0 

Extreme  type  1 

41 1  m/s 

0.10 

[5] 

Geomaterial  bulk  modulus,  K 

Lognormal 

9520  Mpa 

0.15 

[5] 

Density  of  geomaterial,  pg 

Lognormal 

1970  kg/m3 

0.15 

[5] 

Slope  of  linear  fit,  a 

Normal 

0.80 

0.10 

[5] 

Depth  of  target  in  geomaterial,  dg 

Normal 

1.0m 

0.05 

Assumed 

Non-dimensional  parameter,  S 

Normal 

11.10 

0.15 

[2] 

Density  of  concrete,  pc 

Lognormal 

2000  kg/m3 

0.10 

[2] 

Concrete  compressive  strength,  fc 

Lognormal 

40.0  Mpa 

0.10 

[2] 

Thickness  of  concrete  target,  dc 

Normal 

0.50m 

0.05 

Assumed 

Elastic  modulus  of  missile,  E 

Normal 

100  GPa 

0.10 

Depleted  uranium 

1450 
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configuration  is  considered  to  be  nominal  and  a  certain  coefficient  of  variation  (COV)  is  assigned 
for  all  of  the  uncertain  parameters.  This  COV  is  the  measure  of  uncertainty;  and,  the  higher  this 
value,  the  higher  is  the  variation  in  the  parameter. 

Based  on  the  uncertain  variables  defined  in  the  table,  the  reliability  of  the  mission  (penetration 
into  the  concrete  target)  is  determined  to  be  42.37%.  This  would  be  the  final  reliability  of  the 
mission  if  the  missile  was  a  rigid  block.  However,  since  the  missile  is  susceptible  to  buckling  the 
reliability  of  the  missile  is  determined  using  the  second  constraint  discussed  above.  This  value  is 
determined  as  96.38%,  which  means  that  there  is  relatively  low  probability  of  the  missile  buckling 
during  the  mission.  However,  the  system  reliability  is  of  interest  in  this  work  and  this  can  be 
obtained  using  the  above-discussed  methodology.  The  entire  domain  represented  by  uncertain 
variables  is  modeled  using  eight  response  surfaces  that  were  constructed  using  the  Latin 
hypercube  sampling  method  with  60  data  points.  Therefore,  480  data  points  were  required  to 
model  the  response  surfaces  that  were  used  in  the  analysis. 

These  response  surfaces  were  used  to  determine  the  range  of  integrals  in  Eq.  (3)  and  FFT  is  used 
to  evaluate  the  integrals.  By  using  this  method,  the  system  reliability  is  determined  to  be  40.46%. 
There  is  a  5%  drop  in  the  reliability  of  the  system  by  just  considering  the  buckling  constraint. 
Therefore,  this  clearly  indicates  that  the  weapon  system  design  is  a  highly  probabilistic  event  and 
deterministic  techniques  do  not  produce  optimum  configurations. 

In  Fig.  6,  we  can  see  the  PDF  of  the  depth  of  penetration,  which  is  clearly  skewed  to  the  left  of 
2.0.  Due  to  the  shape  of  the  distribution,  a  small  change  in  the  limit  of  the  penetration  depth 
would  move  a  large  area  under  the  PDF  into  the  failure  or  safe  direction,  thereby  making  it  very 
sensitive.  Therefore,  one  can  determine  that  the  PDF  is  most  sensitive  to  the  depths  of  the 
geomaterial  and  concrete,  among  other  things. 

There  are  numerous  parameters  that  can  be  controlled  to  determine  the  configuration  that 
maximizes  the  mission  success.  In  order  to  perform  design  optimization,  information  about  the 
sensitivity  of  mission  success  to  these  parameters  is  important.  The  following  section  discusses  the 
sensitivities  of  various  parameters  to  the  mission  success. 


2500- 
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Fig.  6.  PDF  of  depth  of  penetration. 
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4.  Sensitivity  of  mission  success 

In  order  to  determine  the  parameters  that  are  most  sensitive  to  the  mission  success,  the 
sensitivity  of  the  COY  of  each  of  these  variables  is  examined.  Only  the  parameters  that  result  in  at 
least  5%  variation  in  the  mission  success  rate  are  discussed  below. 


•  Mission  success  vs.  missile  shank  radius : 

Fig.  7  shows  the  relationship  between  the  probability  of  mission  success  and  the  missile  shank 
radius  COV.  Clearly,  as  the  COV  increases  the  probability  of  success  decreases  (about  22%) 
due  to  the  fact  that  there  are  more  realizations  of  a  lower  shank  radius,  which  would  not 
penetrate  into  the  target.  The  lower  shank  radius  results  in  more  failures  due  to  the  violation  of 
the  buckling  constraint.  In  this  particular  case  the  buckling  constraint  would  drive  the  system 
reliability  until  a  point  where  the  shank  radius  is  large  enough  and  that  more  force  is  required 
to  penetrate  than  what  is  available  from  the  impact.  Therefore,  tight  control  of  the  missile 
shank  radius  is  required  in  order  to  maintain  a  high  rate  of  success. 

•  Mission  success  vs.  missile  elastic  modulus : 

Fig.  8  shows  the  rate  of  decrease  of  the  probability  of  mission  success  with  respect  to  the 
variation  in  the  elastic  modulus.  Clearly,  the  8%  decrease  in  this  case  can  also  be  attributed  to 
the  increase  in  the  number  of  missiles  with  lower  elastic  modulus,  which  causes  them  to  buckle 


Fig.  8.  Mission  success  vs.  missile  elastic  modulus. 
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before  the  mission  is  complete.  The  bump  in  the  plot  is  due  to  the  convergence  of  the  system 
probability  solution  at  certain  values  of  COV.  Even  though  modulus  of  elasticity  was 
considered  as  a  random  variable  military  applications  can  dictate  a  tighter  tolerance  due  to 
available  resources.  Therefore,  this  8%  decrease  can  easily  be  avoided  by  specifying  required 
tolerances  on  a  modified  nominal  value.  Since  the  overall  trend  can  be  clearly  seen  from  the 
figure,  more  analyses  were  not  performed  to  recreate  the  curve. 

•  Mission  success  vs.  missile  mass : 

Variation  in  the  missile  mass  resulted  in  about  6.5%  of  variation  in  the  probability  of  success. 
This  did  not  prove  to  be  that  significant  in  the  design  process,  considering  the  variations  due  to 
other  parameters.  Therefore,  any  changes  in  the  mass  of  the  missile  due  to  the  change  in 
contents  would  not  affect  the  mission  drastically.  However,  it  is  still  a  parameter  that  needs  to 
be  monitored  due  to  its  tendency  to  lower  the  mission  success  rate.  The  higher  mass  would 
result  in  larger  kinetic  energy  of  the  missile  as  it  hits  the  ground,  thereby  increasing  the  final 
impact  velocity.  This  impact  velocity  is  responsible  for  the  impact  force  from  conservation  of 
momentum  principle.  Therefore,  when  the  mass  was  increased  the  impact  force  increased  and 
the  weapon  buckled,  leading  to  a  failed  mission  (Fig.  9). 

•  Mission  success  vs.  velocity  of  the  missile : 

It  can  be  seen  from  Fig.  10  and  other  results  that  this  is  the  parameter  that  affects  the  mission 
success  the  most.  There  is  about  a  38%  increase  in  the  rate  of  success  as  the  variation  of  velocity 


Fig.  9.  Mission  success  vs.  mass  of  the  missile. 


Fig.  10.  Mission  success  vs.  velocity  of  the  missile. 
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increases  from  a  low  value  to  high  value.  This  is  due  to  the  fact  that  the  velocity  is  defined  using 
an  extreme  value  distribution  that  had  a  skewness  coefficient  (third  moment)  of  1.14,  which 
means  that  the  data  is  skewed  to  the  right  of  the  mean.  Therefore,  as  the  variation  increases 
there  are  more  and  more  realizations  of  velocity  that  are  on  the  higher  end,  which  leads  to  a 
deeper  penetration.  This  can  be  clearly  seen  in  Fig.  11,  where  PDFs  for  increasing  variation  are 
shown.  The  PDF  that  is  spread  the  most  represents  a  20%  COV,  the  tall  PDF  represents  a  5% 
COV,  and  the  third  PDF  representing  10%  COV.  From  these  PDFs  we  can  see  that  the  higher 
velocities  are  more  probable  when  the  COV  increases,  which  results  in  improved  mission 
success. 

•  Mission  success  vs.  non-dimensional  parameter  “S’”: 

In  this  case,  the  affect  of  a  non-dimensional  parameter  that  is  used  in  the  crater  depth  equation 
is  examined  and  there  is  a  16%  increase  in  the  rate  of  success  as  the  variance  increased. 
Therefore,  this  parameter  also  plays  an  important  role  in  the  mission  success.  This  type  of 
parameter  is  usually  considered  as  an  uncontrollable  parameter  and  its  variation  has  to  be 
determined  with  more  confidence  in  order  to  effectively  utilize  other  parameters  that  are 
controllable.  This  is  a  parameter  that  is  used  as  a  fudge  factor  to  match  the  analytical  results 
with  the  experimental  data.  Therefore,  proper  care  must  be  given  to  determining  this 
parameter,  based  on  its  sensitivity  to  the  mission  success  (Fig.  12). 


Fig.  11.  PDF  of  missile  velocity  for  different  COV. 


Fig.  12.  Mission  success  vs.  non-dimensional  parameter  “*S”. 


1454 


R.C.  Penmetsa  /  International  Journal  of  Mechanical  Sciences  47  (2005)  1442-1454 


5.  Summary 

Based  on  the  analysis  performed,  it  can  be  determined  that  the  overall  mission  success  rate  is 
pretty  low  for  the  configuration  that  was  selected.  However,  parameters  that  drastically  affect  the 
probability  of  mission  success  have  been  identified,  and  modifying  their  mean  values  will  result  in 
a  higher  success.  The  missile  velocity  was  determined  to  have  maximum  positive  sensitivity  to  the 
mission  success  probability  and  the  missile  shank  radius  had  the  maximum  negative  sensitivity. 
Therefore,  a  design  that  has  a  tighter  control  on  the  shank  radius  and  allows  for  higher  velocities 
would  greatly  improve  the  mission  success  probability.  Based  on  the  approach  followed  in  this 
paper,  more  constraints  can  be  added  to  the  study  and  with  the  use  of  efficient  system  reliability 
methods  presented  here;  one  can  determine  the  probability  of  success  of  an  event.  Usually,  a 
mission  involves  multiple  events  that  need  to  occur  in  a  sequence,  therefore,  a  realistic  study 
would  be  one  in  which  all  of  these  events  are  modeled  and  simulated. 
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